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In the course of the XMM-Newton survey of the Small Magellanic Cloud (SMC), 
two new bright X-ray sources were discovered exhibiting the spectral characteris- 
tics of High Mass X-ray Binaries - but revealing only weak evidence for pulsations 
in just one of the objects(at 153s in XMMU J010743.1-715953). The accurate X- 
ray source locations permit the identification of these X-ray source with Be stars, 
thereby strongly suggesting these systems are new Be/X-ray binaries. From blue 
spectra the proposed classification for XMMU J010633. 1-731543 is B0.5-lVe and for 
XMMU J010743.1-715953 it is B2IV-Ve. 

Key words: stars:neutron - X-rays:binaries 



1 INTRODUCTION AND BACKGROUND 

The Be/X-ray systems represent the largest sub-class of all 
High Mass X-ray Binaries (HMXB). A survey of the litera- 
ture reveals that of the ~240 HMXBs known in our Galaxy 
and the Magellanic Clouds (Liu et al., 2005, 2006), >50% 
fall within this class of binary. In fact, in recent years it has 
emerged that there is a substantial population of HMXBs in 
the SMC comparable in number to the Galactic population. 
Though unlike the Galactic population, all except one of the 
SMC HMXBs are Be star systems. In these systems the orbit 
of the Be star and the compact object, presumably a neu- 
tron star, is generally wide and eccentric. X-ray outbursts 
are normally associated with the passage of the neutron star 
close to the circumstellar disk (Okazaki & Negueruela 2001), 
and generally are classified as Types I or II (Stella, White 
& Rosner, 1986). The Type I outbursts occur periodically 



* Based on observations with XMM-Newton, an ESA Science 
Mission with instruments and contributions directly funded by 
ESA Member states and the USA (NASA) 



at the time of the periastron passage of the neutron star, 
whereas Type II outbursts are much more extensive and oc- 
cur when the circumstellar material expands to fill most, or 
all of the orbit. General reviews of such HMXB systems may 
be found in Negueruela (1998), Corbet et al. (2009) and Coe 
et al. (2000, 2009). 

One of the aims of the XMM-Newton large program 
SMC survey (Haberl & Pietsch 2008a) is the ongoing study 
of the Be/X-ray binary population of the SMC, which can be 
used as a star formation tracer for ~5Q Myr old populations 
(Antoniou et al. 2010). Two previous systems have recently 
been identified from this survey: XMMU J005011. 2- 730026 
= SXP 214 (Coe et al., 2011) and XMMU J004814.0- 732204 
= SXP11.87 (Sturm et al 2011). In this paper we present 
the analysis of X-ray and optical data from two newly 
discovered X-ray sources XMMU J010633. 1-731543 and 
XMMU J010743. 1-715953. Both of these sources exhibit 
most of the characteristics of Be/X-ray binary systems, but 
only marginal evidence for X-ray pulsations from the spin- 
ning neutron star in one of the systems. If a weak, or unde- 
tectable pulsed signal is a common phenomenon, then this 
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Figure 1. XMM-Newton EPIC colour image of XMMU J010743. 1-715953 (Left) and XMMU J010633.1-731543 {Right). Red/green/blue 
show logarithmically-spaced intensities in the (0.2 — 1.0)/(1.0 - 2.0)/(2.0 - 4.5) keV band. Circles indicate the extraction regions for 
sources and background events. 



may have significant implications for the number of such sys- 
tems in the Magellanic Clouds and the Milky Way, most of 
which have been primarily identified from their X-ray pul- 
sations. 



analysis, we also created merged time series from all three 
instruments in various energy bands. Light curves, binned 
to have 25 cts bin~^ on average, were created using the SAS 
tasks etimeget and corrected for background, vignetting, 
and point spread function losses by epiclccorr. 



2 X-RAY OBSERVATIONS 

2.1 Observations and data reduction 

The observation setup can be found in Table [T] The two 
X-ray sources were discovered in observations 19 and 27, 
respectively, of the XMM-Newton large program survey of 
the SMC (Haberl et al. 2012, in preparation). XMM-Newton 
has three X-ray telescopes (Aschenbach et al. 2002), one 
equipped with EPIC-pn (Striiuder et al. 2001) and two with 
EPIC-MOS (Turner et al. 2001) CCD detectors in the fo- 
cal plane. Both objects were covered with all three instru- 
ments. We used XMM-Newton SAS ll.O.cQ to process the 
data. Temporal screening was not necessary, since soft pro- 
ton flares were at a low level (i.e. <8 cts ks~^ arcmin"^ for 
EPIC-pn). 

We extracted single and double pixel events from EPIC- 
pn and single to quadruple events from EPIC-MOS, both 
having flag=0. Background events were selected from a 
point source free area on the same CCD for all three instru- 
ments from a circle of radius 30". The source selection re- 
gion was a circle with radius determined by eregionalayse 
to optimise the signal-to-noise ratio. X-ray colour images 
indicating the extraction regions are plotted in Fig. [T] Spec- 
tra and response matrices were created with especget. We 
binned the spectra to a signal-to-noise ratio of 5 for each bin. 
For time series, the photon arrival times were recalculated 
for the solar barycentre. To increase the statistics for timing 

^ Science Analysis Software (SAS), 

|http://xmm.vilspa.esa.es/sas/| 



2.2 X-ray coordinates 

Using emldetect for source detection, we determined posi- 
tions for the X-ray sources in both observations. For obser- 
vation 0601212701 we found ten sources, correlating with 
Chandra X-ray sources from McGowan et al. (2008), which 
we use as astrometric reference. Using this correlation, we 
derived an error weighted average offset of ARA=2.7l"and 
ADec=-0.47". The corrected best-fit position for the 
transient is RA(J2000)=01:06:33.10 and Dec(J2000)=- 
73:15:43.1. The la position uncertainty is 0.57", where we 
assume a systematic uncertainty of 0.5"(Pietsch, Freyberg 
& Haberl 2005) , quadratically added to the statistical uncer- 
tainty. We assign the name XMMU J010633. 1-731543 to this 
source. Analogous for observation 0601211901, we found 25 
sources detected with Chandra, allowing to derive a bore- 
sight correction of ARA=-0.7l"and ADec=1.84", which re- 
sults in a source position of RA(J2000)=01:07:43.13 and 
Dec(J2000)=-71:59:53.5 for XMMU J010743.1-715953. The 
total la position uncertainty is 0.55". 

2.3 Spectral analysis 

The X-ray spectra were modelled with an absorbed power- 
law using xspec (Arnaud 1996) version 12.7.0. The Galactic 
foreground column density was set to NH,gai = 6 x lO'^" 
cm~'^ (Dickey & Lockman 1990) with solar abundances ac- 
cording to Wilms, Allen & McCray (2000). An additional 
column density Nh.shic, having SMC-typical abundances set 
to 0.2 for elements heavier than helium (Russell & Dopita, 
1992), was a free parameter in the fit and accounts for 
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Table 1. XMM-Newton observations. 



ObsID 


Revolution 


Date 


Time 
(UT) 


Instrument 


Filter 


Offax(°) 

['] 


Net Exp. 

[ks] 


Net Cts.C') 


["] 


0601211901 


1827 


2009-11-30 


15:08-23:37 


EPIC-pn 


thin 


11.6 


25.7 


817 


20 


XMMUJ010743.1-715953 






14:45-23:37 


EPIC-MOSl 


medium 


10.9 


31.1 


405 


25 








14:45-23:37 


EPIC-MOS2 


medium 


11.6 


31.1 


400 


25 


0601212701 


1840 


2009-12-26 


07:47-17:39 


EPIC-pn 


thin 


7.5 


30.2 


529 


18 


XMMUJ010633. 1-731543 






07:24-17:39 


EPIC-MOSl 


medium 


8.5 


36.1 


210 


20 








07:24-17:39 


EPIC-MOS2 


medium 


8.5 


36.1 


228 


21 



Off axis angle of the sources. ^''^ Net counts in the (0.2 — 10.0) keV band. ''^^ Radius of the source extraction region. 



photo-electric absorption by the SMC interstellar medium 
and source intrinsic absorption. For each source, the three 
EPIC spectra were fitted simultaneously. Constant factors 
(Cmosi and Cmos2) were allowed to vary to account for in- 
strumental differences, relative to EPIC-pn (Cpn ~ 1). The 
best-fit parameters are listed in Table (2] Uncertainties and 
hmits are given for 90% confidence. The spectra together 
with the best-fit model are plotted in Fig. [21 



2.4 Search for pulsations 

In the case of XMMU J010743. 1-715953, we found an in- 
dication for X-ray pulsations at 6.49 x 10""^ Hz in a Fast 
Fourier Transformation (EFT), as shown in Fig. O By us- 
ing a Bayesian periodic signal detection method (Gregory 
& Loredo 1996), we determined the pulse period with la 
uncertainty to (153.99±0.13) s. The background subtracted 
light curve, convolved with this period, is shown in Fig. [4] A 
test resulted in a maximum of x'^ = 45.7. The significance 
of the signal was estimated to be 2. Oct with a Rayleigh Zf 
test for one harmonic (Haberl & Zavlin 2002, Bucheri et al., 
1983). This assumes a number of 93000 independent trials in 
a frequency range of 10"* to 3 Hz and a maximum Zi of 29. 
For the statistical tests, see Fig. (5] We estimate the pulsed 
fraction to (18.7±4.5)% by assuming a sinusoidal pulse pro- 
file. Further confirmation of this pulse period is necessary, 
demanding higher statistics. 

In the case of XMMU J010633. 1-731543, we could not 
find any indication for pulsations in EFTs. To estimate up- 
per limits for the pulsed fraction, we created pulse profiles in 
the (0.2-10.0) keV band of the merged instrument data. We 
folded the Ught curve at the highest peaks in the according 
EFT at 9.7 x 10"*, 4.0 x 10"^ 1.2 x 10"^ 8.5 x 10"^ and 
3.1 X 10"^ Hz. Analogous to above, we fitted a sine function 
and determined pulsed fractions of 3%, 21%, 23%, 19%, and 
25% with typical uncertainties of ±5% taking into account 
background and using 90% confidence intervals. Therefore 
we estimate the pulsed fraction of XMMU J010633. 1-731543 
to be below ~30% for periods between 0.2 and 10*s. 



2.5 Short term X-ray variabiUty 

The background-corrected instrument-merged hght curves 
are shown in Fig. (6] Evident is the variability of source 
XMMU J010633. 1-731543, but also the light curve of 
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Figure 3. Power density spectrum of source 

XMMU J010743.1-715953 for the merged time series in the 
(0.2-10.0) keV band. 



XMMU J010743. 1-715953 contains some indication of flare- 
like features. A test for variability of the binned fight 
curve resulted in x'^/dof — 101/63 corresponding to a prob- 
ability of constancy of 0.16% for XMMU J010743. 1-715953 
and x'^/dof = 274/38 corresponding to a probability of 



2 X 10" 



for XMMU J010633. 1-731543. 



2.6 Long term X-ray variability 

Source XMMU J010743. 1-715953 was observed 3 times with 
Swift on 2009-12-04, 2011-05-17, and 2011-07-02. We ex- 
tracted spectra within a source region of 30" radius and 
from a background region with 200" radius. Due to short 
exposures of 3.5, 4.3, and 2.6 ks the spectra contain only 33, 
10, and 22 net counts. Using C statistics, we fitted the same 
spectral model as above, only allowing the fiux (an overall 
normalisation) to vary. We derived fluxes of 1.3lQ;g, 3.ll]^;g, 



and 1.2l:ll x 10 



'^^ erg cm ^ s ^ 



consistent (within the 
errors) with the flux given in Table 2 for this source. We 
note, that 91%, 60% and 82% of the net counts are detected 
above 2.0 keV, pointing to a generally hard spectrum of the 
source as seen during the XMM-Newton observation. 

There is no ROSAT source detected within 40" (Haberl 
et al., 2000). In a 20.8 ks ROSAT PSPC observation (Ob- 
sID: rp500250n00) from 1993-10-13 to 1993-10-28 the source 
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Figure 2. EPIC-pn (black), EPIC-MOSl (red), and EPIC-MOS2 (green) spectra of XMMU J010743. 1-715953 (left) and 
XMMU J010633. 1-731543 (right). Top panels show the spectra, together with the best-fit model, lower panels show the residua. 



Table 2. Spectral fitting results. 



source 


H,gal 


H,smc 


r(6) 


Cmosi 


CmOS2 




Lid) 


dof 


^rcd 


XMMU J010743.1-715953 
XMMU J010633.1-731543 


0.6 
0.6 


< 0.6 


0.96l«:i? 

+0.09 
'J- '^-0.09 


-I -17+0.13 
J^-l'-0.12 

i-i'J-o.ie 


-I 10+0.13 

1 in+0.17 
l-l'-'-0.16 


9 O+0.3 
^••^-0.2 


5.65 
0.99 


48 
24 


0.76 
1.42 



Column density of photo-electric absorption by the Galaxy (fixed) and the SMC (free) in 10^^ cm^^. '''' Photon index of the power 
law. t'^' Detected flux in the (0.2-10.0) keV band in 10~13 erg cm"^ s~^. Unabsorbed luminosity in the (0.2-10.0) keV band in 10^^ 

erg s~l, assuming a distance of 60 kpc. 



position was observed ^--^15' offaxis. We derive a 3cr upper 
limit of 32 net cts from a circle of 60" radius around the 
source and a background region of 90". This translates into 
a source flux of 1.1 x 10"^^ erg cm"^ s"^ in the (0.2 - 10.0) 
keV band, by assuming the same spectral model as above 
and using PIMlVlfl Therefore the ROSAT non-detection is 
also consistent with a persistent source. 

The position of source XMMU J010633. 1-731543 was in 
the FoV of one Chandra/ACIS-I observation of the SMC 
wing survey (PI: Coe, ObsID: 5496, MJD 53797). Using the 
CIAO Version 4.C0tool aprates, we derived a 3(t upper limit 
of 1.44 X 10"^ cts s"^ in the (0.2 - 10.0) keV band. Assum- 
ing a time-independent spectrum and using PIMMS, this 
translates into a flux of 3.33 x 10~^^ erg cm^^ s^^. Interest- 
ingly, the Chandra non-detection falls at a time of optical 
brightness (see Section 3.1.1. below). A significant source 
variability by a factor of at least 6.9 proves the transient 
behaviour of this system. 

Similarily, for XMMU J010633. 1-731543, no ROSAT 
match was found within 40" and no upper limit deeper than 
the Chandra value could be derived. 



2 Portable Interactive Multi-Mission Simulator (PIMMS), 
|http: / /heasarc . nasa.gov/Tbors/w3pimms. html 

Chandra Interactive Analysis of Observations (CIAO), 
I http: / / cxc . harvard, edu /ciao/ 1 



3 THE OPTICAL COUNTERPARTS 
3.1 XMMU J010633. 1-731543 

Figure [7| shows the position of the XMM-Newton source 
located upon the OGLE I band image of the region. Un- 
like the image obtained from the DSSII project, this image 
clearly resolves the two candidates labeled as Objects 14 
& 76 (following the fuller OGLE III naming convention of 
SMClll.3.14 and SMClll.3.76). 

3.1.1 Optical & IR photometry 

Figure [8] shows the OGLE-III photometry for the proba- 
ble (Object 14) and possible (Object 76) candidate sources 
to XMMU J010633. 1-731543 is illustrated in Figure [3 It is 
immediately apparent that the two objects exhibit very dif- 
ferent styles of behaviour, and also that one is always much 
fainter than the other. The behaviour of Object 14 is charac- 
teristic of the long-term behaviour patterns of Be stars (see 
Mennickent et al 2002) unlike that of Object 76. Both of 
these light curves were subjected to a Lomb-Scargle period 
analysis. 

The light curve of Object 14 is extremely complex. It 
displays a Type-l/Type-2 light curve, according to the clas- 
sification scheme of Mennickent et al. (2002), i.e. it shows 
both outbursts and high/low states. Such composite light 
curves have only been observed for 2.6% of the stars studied 
by those authors. Over a period of 3000 days the ma:ximum 
amplitude is ~0.7mag in /, which is high compared to most 



Two new X-ray sources in the SMC 5 



15 

Time (ks) 




15 20 
Time (ks) 



Figure 6. Merged corrected light curves from EPIC (0.2 - 10 keV) for XMMUJ010743. 1-715953 (left, MJD 55165) and 
XMMUJ010633.1- 731543 (right, MJD 55191), binned with 470 s and 910 s. 
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Figure 4. X-ray pulse profile of XMMU J010743.1-715953 in var- 
ious energy bands from the merged EPIC time scries. The pulse 
profiles arc background subtracted and normalized to the average 
net count rate of (0.8, 2.8, 1.8, and 5.3) xlQ-^ cts s"!. 
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Figure 5. Top: test for constancy of the light curve of 
XMMU J010743.1-715953 for trial frequencies between 6 and 7 
mHz. Middle: Frequency dependence of the Bayesian odds ratio. 
Bottom: Rayleigh test with one involved harmonic. 

Be stars studied in the SMC. The colour and magnitude 
ranges for Typc-l/Type-2 objects, according to Mennickont 
et al. (2002) are 15.44 <V < 17.01 -0.19 < B-V < -0.00 
and —0.19 < V — I < 0.23. So, the most probable optical 
counterpart of XMMU J010633. 1-731543 is brighter in V 
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Figure 7. OGLE I-band finding ciiart for 
XMMUJ010633. 1-731543. Tlie circle indicates the XMM- 
Newton positional uncertainty. The two candidate OGLE objects 
are indicated as numbers 14 and 76. 
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Figure 8. OGLE I-band photometry for the two candidate 
OGLE objects for XMMUJ010633. 1-731543 shown in Figure[7] 



(see Table [2} than most objects of similar type tabulated by 
Mennickent et al. 

Many difTerent attempts to detrend the data were used 
but none revealed any evidence for coherent behaviour in- 
dicative of either binary motion nor Non Radial Pulsations 
(see Bird, Coe & McBride 2011 for a review of such charac- 
teristics in Be star companions to HMXBs). There is some 
evidence in Years 5 & 6 of different behaviour, but, again, 
no periodic signal was seen. 

On the other hand, the light curve of Object 76 is very 
stable and did reveal some evidence of power at a period 
of 15.97d. Folding the data at this period revealed a weak 
modulation with a possible sinusoidal structure. This kind 
of modulation is often associated with pulsations in stars, 
and the true pulse period may not be 15.97d, but rather 
this may be a harmonic of the true period beating with the 
daily sampling period. 

Thus on both the XMM positional location of 
XMMUJ010633. 1-731543 and on the evidence of the vari- 
ability it seems most Ukely that Object 14 is the correct 
counterpart to the X-ray source. This object is listed as 
number 64194 in the catalogue of Massey (2002). In that 
catalogue the optical magnitudes are listed as B=14.61 and 
V=14.60 but there may be some danger of contamination 
from Object 76 if the two objects were not fully resolved. 
However, from our OGLE III data we have an average V- 
I=-0.11 which, when combined with the Massey V value 
would imply 1=14.71. As can be seen from Figure [8] this 
corresponds to the brightest I band state seen by OGLE III. 
The Massey data were collected on JD2451187 which, un- 
fortunately, pre-dates OGLE III and this object is not in the 
part of the SMC covered by OGLE II. So it is not possible 
to verify that the measurements presented in Massey (2002) 
do not represent some combination of the signals from both 
Objects 14 & 76. However, the colours of (B-V)=0.01 and 



Table 3. I, J, H &; K photometric measurements of the Object 
14 candidate for XMMUJ010633. 1-731543 taken on MJD 52869. 
The B & V magnitudes were taken on MJD 51187 (Massey 2002). 



Waveband 


Magnitude 


Error on Magnitude 


B 


14.61 


0.01 


V 


14.60 


0.01 


I 


15.33 


0.01 


J 


15.52 


0.02 


H 


15.46 


0.02 


K 


15.38 


0.04 



V-I=-0.11 would be typical for the Be stars seen in HMXBs 
in the SMC (Coe et al, 2005). 

Object 14 also falls in the IR Sirius catalogue of the 
SMC (Kato et al, 2007). From there it is possible to obtain 
the IR magnitudes obtained on JD 2452869 which give an 
IR colour of J-K=0.14 - again consistent with the IR colours 
of HMXBs (Coe et al, 2005). AU these optical and IR mag- 
nitudes are tabulated in Table [3] 

The Massey (2002) counterpart, object number 64194, 
with RA=01:06:32.9, Dec=- 73:15:42. 2 (2000), has V = 
14.6 ± 0.02, B ~V = 0.01 ± 0.02, U ~ B = -0.92 ± 0.01, 

V - R = 0.13 ± 0.04. In the MCPS catalogue of Zaritsky 
et al. (2002), the closest counterpart has RA=01:06:32.98, 
Dec=-73: 15:42.5 (2000), which coincides with the Massey 
object. This star has the following magnitudes and colors in 
the MCPS catalogue: U = 13.78±0.034, B = 14.681±0.035, 

V = 15.055 ± 0.037, / = 15.072 ± 0.038 and from 2-MASS, 
J = 15.151 ± 0.059, K' = 14.698 ± 0.10 (as given in the 
MCPS catalogue). Thus the reddening free Q-parameter is 
QuBV = {U - B) - 0.72(_B - V) = -0.93, correspond- 
ing to a Bl star using the Massey colors, and Qubv = 
{U - B)- 0.72(5 -V) = -0.63 using the MCPS magni- 
tudes, corresponding to a B3-4 star. All these spectral types 
are consistent with a Be star. 
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In a 5 arcmin region around this object, the average 
reddening given by Zaritsky et al. (2002) is < Ay >~ 0.26± 
0.28 i.e. E{V - I) = 0.104 ± 0.112. On the other hand, 
from the Haschke et al. (2011) reddening maps, E{V — I) — 
0.052 ± 0.015 (from 5 determinations in regions centered 
within 5 arcmin of the object) and Ay = 3.1 x 0.052 — 0.16. 



3.1.2 Optical spectroscopy 

Spectroscopic observations of both sources discussed in this 
paper were made using telescopes both at the South African 
Astronomical Observatory (SAAO) on 2011 Sep 22, and at 
the European Southern Observatory (ESO) on 2011 Dec 
10. At SAAO the 1.9m telescope of the South African As- 
tronomical Observatory (SAAO) was used solely for red 
spectra. A 1200 lines per mm reflection grating blazed at 
68OOA was used with the SITe CCD which is effectively 
2661798 pixels in size, creating a wavelength coverage of 
62OOA to 6900A . The intrinsic resolution in this mode was 
O.42A /pixel. Data were reduced using Figaro Q. 

At ESO red and blue spectra were taken with the 
ESO Faint Object Spectrograph (EF0SC2) mounted at the 
Nasmyth B focus of the 3.6m New Technology Telescope 
(NTT), La Silla, Chile. The EFOSC2 detector (CCD#40) 
is a Loral/Lesser, Thinned, AR coated, UV flooded, MPP 
chip with 2048x2048 pixels corresponding to 4.l'x4.l'on 
the sky. The instrument was in longslit mode with a slit 
width of 1.5". Grisms 14 and 20 were used for blue and red 
end spectroscopy respectively. Grism 14 has a a grating of 
600 lines mm~ and a wavelength range of AA3095-5085 A 
producing a dispersion of 1 A pixel" ^. The resulting spectra 
have a spectral resolution of ~ 12 A. Grism 20 is one of the 
two new Volume-Phase Holographic grisms recently added 
to EFOSC2. It has 1070 hnes pixel"^ but a smaller wave- 
length range, from 6047-7147 A. This results m a superior 
dispersion of 0.55 A pixeP^ and produces a spectral reso- 
lution for our red end spectra of ~ 6 A. Filter OG530 was 
used to block second order effects. The data were reduced 
using the standard packages available in the Image Reduc- 
tion and Analysis Facility (IRAF). Wavelength calibration 
was implemented using comparison spectra of Helium and 
Argon lamps taken through out the observing run with the 
same instrument configuration. The spectra were normal- 
ized to remove the continuum and a redshift correction was 
applied corresponding to the recession velocity of the SMC. 

The red spectra are shown in Figure [5] and indicate 
only a small decline in the Hq emission over few months 
from an Equivalent Width of -14±lA(SAAO 22 Sep 2011) 
to -10.0±0.2A(ESO 10 Dec 2011). In both cases the line 
profile is consistent with a single peak suggesting that we 
are viewing the circumstellar disk in this system face-on. 

The blue spectra of the OGLE 14 counterpart to the 
XMM object is shown in Figure [10] and it can be used to 
establish the spectral classification of this object. OB stars 
in our own galaxy are classified using the ratio of certain 
metal and helium lines (Walborn & Fitzpatrick 1990) based 
on the Morgan-Keenan (MK; Morgan et al. 1943)) system. 
However, this is unsuitable in lower metallicity environments 
as the metal lines are either much weaker or not present. As 
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Figure 9. Red end spectra of the optical counterpart to 
XMMU J010633.1-731543. The measured equivalent widths of the 
Ha line are -14±lA(22 Sep 2011, upper curve) and -10.0±0.2A(10 
Dec 2011, lower curve). 



such, the optical spectra of XMMUJ010632. 4- 731542 and 
XMMUJ010743. 3-715955 were classified using the method 
developed by Lennon (1997) for B-type stars in the SMC and 
implemented for the SMC, LMC and Galaxy by Evans et al 
(2004, 2006). This system is normalized to the MK system 
such that stars in both systems show the same trends in their 
line strengths. The luminosity classification method from 
Walborn & Fitzpatrick (1990) was assumed in this work. 

The optical spectrum of 

XMMU J010633. 1-731543 shows evidence of Siiv lines 
at A4088 and A4116 implying the star cannot be later than 
type Bl. The Hell lines are difficult to distinguish above the 
noise level of the spectrum, however the presence of a line 
at A4686 and the apparent lack of any line at A4541 would 
suggest a spectral classification of BO. 5, as such we classify 
XMMU J010633. 1-731543 as a B0.5-le star. The luminosity 
class was determined using the ratios of S HI A4553 to H I 
A4387 and He I A4143 to He I A4121. The latter increases 
dramatically with increasing luminosity class (from I-V) 
whereas the former decreases with luminosity class. The 
relative strengths of these lines suggests a luminosity 
class of V for XMMU J010633. 1-731543. So the proposed 
classification for XMMU J010633. 1-731543 is B0.5-lVe. 

This spectral classification can be checked against the 
photometry. If we assume that the reddening and bright- 
ening is mainly due to the contribution from the cooler 
equatorial disk, then we can take the color correspond- 
ing to the faintest magnitude as corresponding to the star, 
reddened basically only by interstellar extinction. For the 
faintest V mag recorded, V — I — —0.124, giving (V — I)o = 
{V - I)- E{V - I) = -0.124 - 0.052 = -0.176 which is 
consistent with a B6 MS star (or using the Zaritsky redden- 
ing, (F-7)o=-0. 124-0.104^-0. 23 consistent with a B4-5 MS 
star). The V magnitude corresponding to the bluest (low 
state) color is V = 15.19, i.e. Vo = 15.19 - 0.16 = 15.03, 
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Figure 10. Blue end spectra of the optical counterparts to XMMU J010633. 1-731543 and XMMU J010743.1-715953. 



and for the distance modulus of the SMC of 18.90 ± 0.18 
(Kapakos et al. 2011) Mv = -3.87 consistent with a Bl MS 
star, where we have used the Haschke et al. reddening value. 
For the Zaritsky et al. reddening value, the corresponding 
value is Mv = —3.97 more consistent with a BO MS star. 

So both the Q-parameter and the above reasoning lead 
the photometric data to suggest a B0-B5 main sequence star, 
consistent with the spectral classification. 



3.2 XMMU J010743. 1-715953 

Figure [TT] shows the position of the XMM source located 
upon the DSSII R band image of the region. One object 
clearly lies in the XMM uncertainty region and is indi- 
cated in this figure - the OGLE III object SMC 113.4.6850. 
This object is also known as [MA93] 1640 (Meyssonnier 
& Azzopardi, 1993). They describe the H-alpha emission 
of this object as weak but sharp. The candidate is also 
to be found in the 2MASS catalogue as 2MASS01074325- 
7159539 with J = 16.293 ± 0.101, H = 16.253 ± 0.250 and 
K = 15.839 ± 0.259. It is also recorded in the DENIS cat- 
alogue of Cioni et al. 2000, with I = 16.358 ± 0.07 and 
J = 16.591 ±0.24. 




Figure 11. DSSII R-band finding chart (5x5 arcmin) for 
XMMU J010743.1-715953. The position of the optical counter- 
part is indicated. 

3.2.1 Optical & IR photometry 

Figure [12] shows the OGLE-III I-band photometry for the 
proposed candidate source to XMMU J010743. 1-715953 il- 
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Figure 12. OGLE I-band photometry for the 
XMMUJ010743. 1-715953 candidate OGLE object SMC 
113.4.6850 shown in Figure [TT] 
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Figure 13. Red end spectra of the optical counterpart to 
XMMU J010743.1-715953. The measured equivalent widths of the 
Ha line are -14±lA(14 Dec 2009), -5±5A(22 Sep 2011) and 
-|-1.5±0.5A(10 Dec 2011). 



Table 4. I, J, H & K photometric measurements of the 
XMMU J010743. 1-715953 candidate - see text for details. 



Waveband 


Magnitude 


Magnitude 




MJD52943 


MJD55179 


I 


16.56±0.01 




J 


16.40±0.02 


16.36±0.05 


H 


16.21±0.02 


16.23±0.09 


K 


16.00±0.06 


16.03±0.22 



lustrated in Figure [TT] It is immediately apparent that there 
are significant brightness changes in the system typical of Be 
star behaviour. 

The candidate for XMMU J010743. 1-715953 also falls 
in the IR catalogue of the SMC (Kato et al, 2007) with the 
name Sirius 01074343-7159539. From there it is possible to 
obtain the IR magnitudes obtained on JD 2452943 (30 Oct 
2003) which give an IR colour of J-K=0.40 - again very con- 
sistent with the IR colours of HMXBs (Coe et al, 2005). It is, 
however, rather faint with its J magnitude being the faintest 
of all know Be/X-ray binary systems in the SMC. The IR 
values were confirmed by further observations with the same 
camera and telescope six years later on JD 2455179 (14 Dec 
2009). In addition to the IR measurements, an OGLE-III 
I-band measurement exists for the same date as the first IR 
data set. All these magnitudes are tabulated in Table [T] 

3.2.2 Optical spectroscopy 

Spectroscopic observations were taken from both SAAO 
(2009 Dec 14 & 2011 Sep 22) and ESQ (2011 Dec 10) as 
described above in Section 3.1.2. 

In comparison to XMMU J010633. 1-731543, the blue 
optical spectrum XMMU J010743. 1-715953 shown in Fig- 
ure [To] appears featureless around the H5 line. There is no 



indication of any He II or Si IV lines suggesting a spec- 
tral classification later than B1.5. The presence of the Si 
III line at A4553 sets a limit of B3 on the spectral clas- 
sification. Stars of spectral type B2-3 are classified based 
on the relative strength of the Si III with respect to the 
Mg II line at A4481. It is difficult to determine if this line 
is even present, however we can rule out the case of a 
stronger Mgll A4481 with respect to Si III excluding a clas- 
sification of B3. We tentatively classify XMMUJ010743.3- 
715955 as a B2e star. As for the previos object, the lumi- 
nosity class was determined using the ratios of S III A4553 
to H I A4387 and He I A4143 to He I A4121. The relative 
strengths of these lines suggests a luminosity class of IV - 
V for XMMU J010743. 1-715953. So the proposed spectral 
classification for XMMU J010743.1-715953is B2IV-Ve. 

The red spectra are shown in Figure [13] It is clear that 
Ha line evolved from a strong, double-peaked emission fea- 
ture to an absorption line showing little or no evidence for 
the presence of a circumstellar disk. This is supported by 
the lack of any H/3 emission in Figure [lO] 



4 DISCUSSION AND CONCLUSIONS 
4.1 XMMU J010633. 1-731543 

The X-ray spectrum of XMMU J010633. 1-731543 shows no 
significant X-ray absorption in addition to the Galactic one, 
suggesting that this system is located on the near side of the 
SMC and also is probably viewed face-on. The total SMC HI 
column density in this direction is 3.95 x 10^^ cm~^ . The lack 
of pulsations in this object can be caused by many possible 
reasons e.g. weak magnetic fields, geometric effects, excep- 
tional long pulse periods, or simply insufficient statistics on 
lower luminosity sources. 

It is very clear from Figure [8] that there are epochs of 



10 M.J. Coe et al. 




wind outflow, and then the slopes shown indicate the grad- 
ual dissipation of the circumstellar disk all with the same 
characteristic timescale. This idea would be supported by 
the colour changes seen at the same time. 

4.2 XMMUJ010743. 1-715953 



15.4 15.2 15 14.8 14.6 

I band magnitude 

Figure 14. Colour changes observed in the OGLE 14 optical 
counterpart to XMMU J010633.1-731543 in the period MJD3300- 
4800. 



dramatic change in the brightness of OGLE SMClll.3.14. 
In particular around MJD 54200-54600 there is a consider- 
able decrease in the I band magnitude in an linear manner 
on this diagram. The slope of the decrease is ~0.5 mag in 
130d. Similar slopes may be seen around MJD 52800-52900 
and MJD 2200. In all cases the rate of change is very similar 
to the large change. To explore this further we used avail- 
able OGLE III V-band data to study any colour changes 
that may be occurring. For this source there are OGLE III 
V-band measurements available for the period MJD53400- 
54800 which includes the largest variation shown in Figure [H] 
Using these data it is possible to construct a magnitude- 
colour diagram and this is presented in Figure [ill Note 
colour information is only available on occasions when the I 
and V band data were taken within 24 hours of each other. 
It is immediately apparent that there are substantial colour 
changes taking place as the source fades and brightens, in 
the sense that the source becomes bluer as it fades (lower 
right hand corner to top left corner of the figure). 

Similar color-magnitude correlations in BeXRBs have 
also been observed in the past, e.g. in McGowan & Charles 
2002 (MNRAS), Clark et al.(1999), and more recently in 
Sarty G.E. et al. (2011). 

Janot-Pacheco et al. (1987) have shown that changes 
in the emission from the circumstellar disk will produce a 
correlation between V magnitude and B-V color, with the 
system becoming redder as it brightens owing to the larger 
contribution from the cooler disk. 

An obvious interpretation for this is that the decline of 
the circumstellar disk inevitably takes away cooler, redder 
light from the overall emission of the system. 

Since there are at least 3 epochs of major declines in 
brightness (and a couple of shorter ones) shown in Figure [8] 
it is obviously a phenomenon characteristic of this source. 
One possible explanation could be that the start of each of 
these declines coincides with the switching off of the stellar 



For source XMMU J010743.1-715953, the total SMC HI col- 
umn density is 6.3 xlO'^^ cm~'^, a factor of ~12.8 smaller 
than the value derived from the X-ray absorption. Therefore 
the high X-ray absorption of XMMU J010743. 1-715953 has 
to be source intrinsic. Also compared to other BeXRB, the 
value is quite high (Haberl, Eger & Pietsch, 2008)and the 
photon distribution measured with Swift is rather hard. One 
possible scenario is that the neutron star is seen in/through 
the circumstellar disk. Certainly the double peaked Ha line 
seen in Figure [131 indicates a high inclination to the system. 

The Swift detection of XMMU J010743. 1-715953 four 
days after the XMM-Newton observation would be in agree- 
ment with a Type I outburst. In contrast to that, the two 
Swift detections in 2011 are separated by 46 days, and have 
a rather low luminosity and are therefore unlikely caused by 
a Type II outburst. The fact that the source is detected at 
3 epochs at a rather low and similar luminosity suggest a 
persistent BeXRB. This would be expected, if the neutron 
star orbit is in the plane of the circumstellar disk of the Be 
star and has a low eccentricity. 

The timescale of the optical variations seen in Figure [T^ 
is of the order of ten years or more. The overall shape of the 
light-curve can probably be classified as a Type-3 light curve, 
according to Mennickent et al (2002). However it must be 
noted that all of the light-curves that they classify as Type-3 
have shorter characteristic timescales of modulation (of the 
order of 1 yr). It is possible (see discussion on Type-3 light- 
curves in their paper) that this very long-term periodic or 
semi-periodic variability is linked to global one-armed oscil- 
lations in the circumstellar envelope (Okazaki et al., 1997). 
The period of the modulation is about 7 years on average 
and it is a weak function of spectral type, increasing as the 
spectral type becomes later (see Fig. 4 in Okazaki et al 1997). 
The associated photometric variations are not strictly peri- 
odic. The problem is that these oscillations are not expected 
theoretically to occur in the low metallicity environment of 
the SMC. However, Hummel et al. (2001) found observa- 
tional evidence contradicting this theoretical prediction. 

In addition to the general changes in fiux, the presence 
of sharp features can be seen, in particular when the source 
is in a fainter state around MJD 4000-4400. If the OGLE 
HI data are detrended and then subject to Lomb-Scargle 
analysis a strong peak emerges in the power spectrum at 
a period of 100.3d. This period is present in both the first 
half and the last half of the data, and not just a feature of 
the later years. Though it is probable that the depth of the 
modulation increases when the source is generally fainter. 
The result of folding the last three years of I-band data at 
the period of 100.3d is shown in the upper panel of Figure lTSl 

Shown in Figure [TBI are the OGLE III I-band data con- 
verted to flux units with the global trend removed. It is ap- 
parent from this flgure that the depth of the dips and spikes 
grows in size throughout the data set. A result with some 
similarities to these data has been presented by Graczyk et 
al, 2011 also from OGLE III data of the Magellanic Clouds. 
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Figure 15. The last three years of detrended OGLE I-band data 
folded at the period of 100.3d and with a ro=MJD3872.0 for the 
optical counterpart to XMMUJ010743. 1-715953. The vertical 
scale is detrended magnitudes. 



Figure 17. OGLE data for XMMU J010743.1-715953 showing 
the colour change pattern. Each ellipse approximately encloses 
data from an annual observing season - each season has a different 
symbol. The total period covered by the above graph (MJD 3326 
- 4953) corresponds to the major change in optical brightness in 
the source shown in Figure [12] 
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relative positions of each year's data set it can be seen that 
the colour of the system follows and anti-clockwise pattern 
on this diagram. Similar colour evolutions seen in other SMC 
Be/X-ray binary systems have been presented in Rajoeli- 
manana et al (2011). Starting in the top right hand corner 
of Figure [17] it can be seen that the source fades in overall 
brightness by about 0.3 magnitude with only small changes 
in colour, but then as the disk re-builds (top left to lower 
right) there is a significant increase in the redness of the 
system. This presumably reflects the growth of the disk as 
a large cooler component of the system. In both cases the 
timescales for the disk loss and subsequent rebuilding is 1-2 
years. 



Figure 16. Detrended OGLE I-band data for 
XMMU J010743. 1-715953 expressed in flux units. 

In their Figure 9 they show an example of a Transient Eclips- 
ing Binary system in the LMC which shows strong evidence 
for variable-depth, eclipse-like features. These authors ex- 
plain the variation in the eclipse features as due to regres- 
sion of the nodes of the orbital plane and compare the result 
to the galactic system V907 Sco (Lacy et al, 1999). Such re- 
gression is attributed to a third body in the system. 

Shown in Figure[T7]is the evolution of the colours of the 
optical counterpart to XMMU J010743. 1-715953 seen during 
the major brightness change shown in Figure [121 From the 
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